Editorial : Sources and propagation of ultra-low frequency waves in planetary magnetospheres

Non peer reviewe

Statistical study of the alteration of the magnetic structure of magnetic clouds in the Earth’s magnetosheath

The magnetosheath plays a central role in the solar wind-magnetospheric coupling. Yet the effects of its crossing on solar wind structures such as magnetic clouds (MCs) are generally overlooked when assessing their geoeffectivity. Using 82 MCs observed simultaneously in the solar wind and the magnetosheath, we carry out the first statistical study of the alteration of their magnetic structure in the magnetosheath. For each event, the bow shock properties are obtained from a magnetosheath model. The comparison between the model results and observations shows that in 80% of cases, the MHD-based model captures well the magnetosheath transition; the other events are discussed separately. We find that just downstream of the bow shock the variation of the magnetic field direction shows a very good anticorrelation (r = −0.91) with the angle between the upstream magnetic field and the shock normal. We then focus on the magnetic field north-south component Bz because of its importance for geoeffectivity. Although the sign of Bz is generally preserved in the magnetosheath, we also find evidence of long-lasting intervals of opposite Bz signs in the solar wind and the magnetosheath during some events, with a |Bz| reversal > 10 nT at the magnetopause. We find that these reversals are due to the draping of the field lines and are associated with predominant upstream By. In those cases, the estimated position of the regions of antiparallel fields along the magnetopause is independent of the sign of the upstream Bz . This may have strong implications in terms of reconnection.Peer reviewe

Non-locality of Earth's quasi-parallel bow shock : injection of thermal protons in a hybrid-Vlasov simulation

We study the interaction of solar wind protons with Earth's quasi-parallel bow shock using a hybrid-Vlasov simulation. We employ the global hybrid model Vlasiator to include effects due to bow shock curvature, tenuous upstream populations, and foreshock waves. We investigate the uncertainty of the position of the quasi-parallel bow shock as a function of several plasma properties and find that regions of non-locality or uncertainty of the shock position form and propagate away from the shock nose. Our results support the notion of upstream structures causing the patchwork reconstruction of the quasi-parallel shock front in a non-uniform manner. We propose a novel method for spacecraft data to be used to analyse this quasi-parallel reformation. We combine our hybrid-Vlasov results with test-particle studies and show that proton energization, which is required for injection, takes place throughout a larger shock transition zone. The energization of particles is found regardless of the instantaneous non-locality of the shock front, in agreement with it taking place over a larger region. Distortion of magnetic fields in front of and at the shock is shown to have a significant effect on proton injection. We additionally show that the density of suprathermal reflected particles upstream of the shock may not be a useful metric for the probability of injection at the shock, as foreshock dynamics and particle trapping appear to have a significant effect on energetic-particle accumulation at a given position in space. Our results have implications for statistical and spacecraft studies of the shock injection problem.Peer reviewe

Hybrid-Vlasov simulation of auroral proton precipitation in the cusps : Comparison of northward and southward interplanetary magnetic field driving

Particle precipitation is a central aspect of space weather, as it strongly couples the magnetosphere and the ionosphere and can be responsible for radio signal disruption at high latitudes. We present the first hybrid-Vlasov simulations of proton precipitation in the polar cusps. We use two runs from the Vlasiator model to compare cusp proton precipitation fluxes during southward and northward interplanetary magnetic field (IMF) driving. The simulations reproduce well-known features of cusp precipitation, such as a reverse dispersion of precipitating proton energies, with proton energies increasing with increasing geomagnetic latitude under northward IMF driving, and a nonreversed dispersion under southward IMF driving. The cusp is also found more polewards in the northward IMF simulation than in the southward IMF simulation. In addition, we find that the bursty precipitation during southward IMF driving is associated with the transit of flux transfer events in the vicinity of the cusp. In the northward IMF simulation, dual lobe reconnection takes place. As a consequence, in addition to the high-latitude precipitation spot associated with the lobe reconnection from the same hemisphere, we observe lower-latitude precipitating protons which originate from the opposite hemisphere's lobe reconnection site. The proton velocity distribution functions along the newly closed dayside magnetic field lines exhibit multiple proton beams travelling parallel and antiparallel to the magnetic field direction, which is consistent with previously reported observations with the Cluster spacecraft. In both runs, clear electromagnetic ion cyclotron waves are generated in the cusps and might further increase the calculated precipitating fluxes by scattering protons to the loss cone in the low-altitude cusp. Global kinetic simulations can improve the understanding of space weather by providing a detailed physical description of the entire near-Earth space and its internal couplings.Peer reviewe

Phase space density analysis of outer radiation belt electron energization and loss during geoeffective and nongeoeffective sheath regions

Coronal mass ejection driven sheath regions are one of the key drivers of drastic outer radiation belt responses. The response can however be significantly different based on the sheath properties and the associated inner magnetospheric wave activity. We performed two case studies on the effects of sheaths on outer belt electrons of various energies using data from the Van Allen Probes. One sheath caused a major geomagnetic disturbance and the other had only a minor impact. We especially investigated the phase space density (PSD) of seed, core, and ultrarelativistic electrons to determine the dominant energization and loss processes taking place during the events. Both sheaths produced substantial variation in the electron fluxes from tens of kiloelectronvolts up to ultrarelativistic energies. The responses were however the opposite: the geoeffective sheath mainly led to enhancement, while the nongeoeffective one caused a depletion throughout most of the outer belt. The case studies highlight that both inward and outward radial transport driven by ultra-low frequency waves played an important role in both electron energization and loss. Additionally, PSD radial profiles revealed a local peak that indicated significant acceleration to core energies by chorus waves during the geoeffective event. The distinct responses and different mechanisms in action during these events were related to the timing of the peaked solar wind dynamic pressure causing magnetopause compression, and the differing levels of substorm activity. The most remarkable changes in the radiation belt system occurred in key sheath sub-regions near the shock and the ejecta leading edge.Peer reviewe

Transmission of an ICME Sheath Into the Earth's Magnetosheath and the Occurrence of Traveling Foreshocks

The transmission of a sheath region driven by an interplanetary coronal mass ejection into the Earth's magnetosheath is studied by investigating in situ magnetic field measurements upstream and downstream of the bow shock during an ICME sheath passage on 15 May 2005. We observe three distinct intervals in the immediate upstream region that included a southward magnetic field component and are traveling foreshocks. These traveling foreshocks were observed in the quasi-parallel bow shock that hosted backstreaming ions and magnetic fluctuations at ultralow frequencies. The intervals constituting traveling foreshocks in the upstream survive transmission to the Earth's magnetosheath, where their magnetic field, and particularly the southward component, was significantly amplified. Our results further suggest that the magnetic field fluctuations embedded in an ICME sheath may survive the transmission if their frequency is below similar to 0.01 Hz. Although one of the identified intervals was coherent, extending across the ICME sheath and being long-lived, predicting ICME sheath magnetic fields that may transmit to the Earth's magnetosheath from the upstream at L1 observations has ambiguity. This can result from the strong spatial variability of the ICME sheath fields in the longitudinal direction, or alternatively from the ICME sheath fields developing substantially within the short time it takes the plasma to propagate from L1 to the bow shock. This study demonstrates the complex interplay ICME sheaths have with the Earth's magnetosphere when passing by the planet.Peer reviewe

Transmission of an ICME Sheath Into the Earth's Magnetosheath and the Occurrence of Traveling Foreshocks

The transmission of a sheath region driven by an interplanetary coronal mass ejection into the Earth's magnetosheath is studied by investigating in situ magnetic field measurements upstream and downstream of the bow shock during an ICME sheath passage on 15 May 2005. We observe three distinct intervals in the immediate upstream region that included a southward magnetic field component and are traveling foreshocks. These traveling foreshocks were observed in the quasi-parallel bow shock that hosted backstreaming ions and magnetic fluctuations at ultralow frequencies. The intervals constituting traveling foreshocks in the upstream survive transmission to the Earth's magnetosheath, where their magnetic field, and particularly the southward component, was significantly amplified. Our results further suggest that the magnetic field fluctuations embedded in an ICME sheath may survive the transmission if their frequency is below similar to 0.01 Hz. Although one of the identified intervals was coherent, extending across the ICME sheath and being long-lived, predicting ICME sheath magnetic fields that may transmit to the Earth's magnetosheath from the upstream at L1 observations has ambiguity. This can result from the strong spatial variability of the ICME sheath fields in the longitudinal direction, or alternatively from the ICME sheath fields developing substantially within the short time it takes the plasma to propagate from L1 to the bow shock. This study demonstrates the complex interplay ICME sheaths have with the Earth's magnetosphere when passing by the planet.Peer reviewe

Outer radiation belt and inner magnetospheric response to sheath regions of coronal mass ejections : a statistical analysis

The energetic electron content in the Van Allen radiation belts surrounding the Earth can vary dramatically at several timescales, and these strong electron fluxes present a hazard for spacecraft traversing the belts. The belt response to solar wind driving is, however, largely unpredictable, and the direct response to specific large-scale heliospheric structures has not been considered previously. We investigate the immediate response of electron fluxes in the outer belt that are driven by sheath regions preceding interplanetary coronal mass ejections and the associated wave activity in the inner magnetosphere. We consider the events recorded from 2012 to 2018 in the Van Allen Probes era to utilise the energy- and radial-distance-resolved electron flux observations of the twin spacecraft mission. We perform a statistical study of the events by using the superposed epoch analysis in which the sheaths are superposed separately from the ejecta and resampled to the same average duration. Our results show that the wave power of ultra-low frequency Pc5 and electromagnetic ion cyclotron waves, as measured by a Geostationary Operational Environmental Satellite (GOES), is higher during the sheath than during the ejecta. However, the level of chorus wave power, as measured by the Van Allen Probes, remains approximately the same due to similar substorm activity during the sheath and ejecta. Electron flux enhancements are common at low energies ( 4). It is distinctive that the depletion extends to lower energies at larger distances. We suggest that this L-shell and energy-dependent depletion results from the magnetopause shadowing that dominates the losses at large distances, while the wave-particle interactions dominate closer to the Earth. We also show that non-geoeffective sheaths cause significant changes in the outer belt electron fluxes.Peer reviewe

On the Importance of Spatial and Velocity Resolution in the Hybrid-Vlasov Modeling of Collisionless Shocks

In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role.In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role.In hybrid-Vlasov plasma modeling, the ion velocity distribution function is propagated using the Vlasov equation while electrons are considered a charge-neutralizing fluid. It is an alternative to particle-in-cell methods, one advantage being the absence of sampling noise in the moments of the distribution. However, the discretization requirements in up to six dimensions (3D position, 3V velocity) make the computational cost of hybrid-Vlasov models higher. This is why hybrid-Vlasov modeling has only recently become more popular and available to model large-scale systems. The hybrid-Vlasov model Vlasiator is the first to have been successfully applied to model the solar-terrestrial interaction. It includes in particular the bow shock and magnetosheath regions, albeit in 2D-3V configurations so far. The purpose of this study is to investigate how Vlasiator parameters affect the modeling of a plasma shock in a 1D-3V simulation. The setup is similar to the Earth's bow shock in previous simulations, so that the present results can be related to existing and future magnetospheric simulations. The parameters investigated are the spatial and velocity resolution, as well as the phase space density threshold, which is the key parameter of the so-called sparse velocity space. The role of the Hall term in Ohm's law is also studied. The evaluation metrics used are the convergence of the final state, the complexity of spatial profiles and ion distributions as well as the position of the shock front. In agreement with previous Vlasiator studies it is not necessary to resolve the ion inertial length and gyroradius in order to obtain kinetic phenomena. While the code remains numerically stable with all combinations of resolutions, it is shown that significantly increasing the resolution in one space but not the other leads to unphysical results. Past a certain level, decreasing the phase space density threshold bears a large computational weight without clear physical improvement in the setup used here. Finally, the inclusion of the Hall term shows only minor effects in this study, mostly because of the 1D configuration and the scales studied, at which the Hall term is not expected to play a major role.Peer reviewe

Asymmetries in the Earth's dayside magnetosheath : results from global hybrid-Vlasov simulations

Bounded by the bow shock and the magnetopause, the magnetosheath forms the interface between solar wind and magnetospheric plasmas and regulates solar wind-magnetosphere coupling. Previous works have revealed pronounced dawn-dusk asymmetries in the magnetosheath properties. The dependence of these asymmetries on the upstream parameters remains however largely unknown. One of the main sources of these asymmetries is the bow shock configuration, which is typically quasi-parallel on the dawn side and quasi-perpendicular on the dusk side of the terrestrial magnetosheath because of the Parker spiral orientation of the interplanetary magnetic field (IMF) at Earth. Most of these previous studies rely on collections of spacecraft measurements associated with a wide range of upstream conditions which are processed in order to obtain average values of the magnetosheath parameters. In this work, we use a different approach and quantify the magnetosheath asymmetries in global hybrid-Vlasov simulations performed with the Vlasiator model. We concentrate on three parameters: the magnetic field strength, the plasma density, and the flow velocity. We find that the Vlasiator model reproduces the polarity of the asymmetries accurately but that their level tends to be higher than in spacecraft measurements, probably because the magnetosheath parameters are obtained from a single set of upstream conditions in the simulation, making the asymmetries more prominent. A set of three runs with different upstream conditions allows us to investigate for the first time how the asymmetries change when the angle between the IMF and the Sun-Earth line is reduced and when the Alfven Mach number decreases. We find that a more radial IMF results in a stronger magnetic field asymmetry and a larger variability of the magnetosheath density. In contrast, a lower Alfven Mach number leads to a reduced magnetic field asymmetry and a decrease in the variability of the magnetosheath density, the latter likely due to weaker foreshock processes. Our results highlight the strong impact of the quasi-parallel shock and its associated foreshock on global magnetosheath properties, in particular on the magnetosheath density, which is extremely sensitive to transient quasi-parallel shock processes, even with the perfectly steady upstream conditions in our simulations. This could explain the large variability of the density asymmetry levels obtained from spacecraft measurements in previous studies.Peer reviewe